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We study the sensitivity of the r-process abundance pattern to neutron capture rates along the 
rare earth region (^4 ~ 150 to A ~ 180). We introduce the concepts of large nuclear flow and flow 
saturation which determine the neutron capture rates that are influential in setting the rare earth 
abundances. We illustrate the value of the two concepts by considering high entropy conditions 
favorable for rare earth peak production and identifying important neutron capture rates among 
the rare earth isotopes. We also show how these rates influence nuclear flow and specific sections of 
the abundance pattern. 
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I. INTRODUCTION 

The rapid neutron capture process or 'r-process' has 
long been known to be an integral component of heavy 
element nucleosynthesis. The onset of the r-process has 
traditionally been characterized by a relatively large neu- 
tron number density (n n > 10 20 cm~ 3 ) and high temper- 
ature (T~ 10GK). These quantities decrease culminating 
in the last stage of the r-process known as (n, 7) <=± (7, n) 
freeze-out [T]. 

During freeze-out nuclides fall out of out of (71,7) 
(7, n) equilibrium and individual neutron capture rates 
are important in determining final abundances [2H5]. 
For a classical, 'hot' freeze-out temperatures around 
T> 1GK are expected. Recently, Wanajo jS] has sug- 
gested that r-process nuclides may participate in a 'cold' 
(n, 7) <=2 (7, n) freeze-out, with temperatures as low as 
(T~ .1GK). Both hot and cold freeze-out scenarios fa- 
vor short timescales for neutron capture which can first 
exceed and then compete with /3-decay. Nucleosynthe- 
sis in these environments progresses along the NZ-plane, 
traversing the nuclear landscape including the rare earth 
region (A ~ 160) far from the valley of beta stability 
where little experimental nuclear data exists [7J. 

The astrophysical location of the r-process is not 
known at this time [IJ [8] [5]. There are several candi- 
date sites where the r-process may occur. Among the 
possible candidates are the high entropy senario of su- 
pernova ejecta (neutrino driven wind) [104414] . supernova 
fallback [15], collapse of O-Ne-Mg cores [TBlUS] . gamma 
ray burst accretion disks |19H21| . and various compact 
object merging scenarios [22H24] . Observational data 
from metal-poor stars [2SH27J favors massive stars which 
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mature on short timescales. A comparison of sites has 
suggested that core-collapse supernova should be favored 
over neutron star mergers for the production of the heav- 
iest r-process elements [35]. However, modern supernova 
calculations typically do not achieve sufficiently neutron 
rich conditions conducive to heavy element nucleosynthe- 
sis [DM]. 

Our understanding of the r-process is also naturally en- 
tangled with our knowledge of the input nuclear physics. 
For example, it is well known that the long /3-decay half- 
lifes found at closed neutron shells are responsible for the 
A = 130 and A — 195 peaks found in the r-process abun- 
dance pattern [52] . see Figure [l] While it is difficult to 
measure the properties of short-lived nuclides far from 
stability there have been many recent advances [33h40j . 
Future radioactive beam facilities are expected to help in 
this endeavor. Due to limited data on nuclides far from 
stability, theoretical extrapolations must be used as input 
for r-process calculations. 

The primary nuclear physics inputs needed to deter- 
mine the nucleosynthetic outcome include neutron cap- 
ture cross sections, separation energies, masses, /3-decay 
rates, and branching ratios. The effects of different 
masses [H] and /3-decay rates [32 03] have been studied 
in detail for some time. Until recently however, neutron 
capture rate cross sections have warranted less consider- 
ation. 

In the past decade studies of neutron capture rates 
have been performed by several groups C3 S3H3Z] ■ Two 
recent studies of particular interest have been performed 
by Beun et al. [2], and Surman et al. [3] who demon- 
strated capture rate effects of individual isotopes. The 
former's study focused on the single neutron capture of 
130 Sn. The effects of elements near the A = 130 peak 
were studied in detail by Surman et al. [3j. In each of 
these two studies the neutron capture rates of single el- 
ements in the A = 130 region lead to global changes in 
the final r-process abundance pattern. 

It has been shown in [5J 48 that neutron capture rates 
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are important for correctly forming the rare earth peak. 
In order to understand observational data e.g. [49], 150]. 
it is crucial to determine which rates are important and 
how these rates influence nuclear flow and ultimately final 
abundances. If in the future improved calculations or 
measurements are included this will increase the efficacy 
of the rare earth peak as a freeze-out diagnostic [5T] . 

In this paper we study neutron capture rates in the rare 
earth region using ideal freeze-out conditions for this sec- 
tor [5J [51] . Modifying individual neutron capture rates 
of nuclides in this region, we identify the most influen- 
tial and in return show their effects on sections of the re- 
process abundance pattern. We also isolate factors which 
influence the magnitude of this neutron capture rate ef- 
fect. 



II. CALCULATIONS 

To study the effects of neutron capture rates on final 
abundances we implement a one-dimensional model of 
the r-process by following the abundance composition of 
a single ejected mass zone. It was shown in [5j 5Tj that 
the rare earth elements are sensitive to the thermody- 
namic evolution of the ejected material. Thus, to sim- 
ulate the qualitative behavior of hydrodynamic outflows 
we parameterize our trajectories via two analytic proce- 
dures. For the hot freeze-out evolution we parameterize 
the density with the same functional form as [53) . 



P (t) = Pl CX P {-t/T) + P2 ^^- t ) (1) 

where p\ + pi is the density at time t — 0, 3r = Td yn and 
A chosen so that the two terms on the right hand side 
are equivalent at a time t = r. The first term controls 
the trajectory at early times, during which the neutron- 
to-seed ratio is set and the second term controls the late 
time behavior, during which the rare earth peak forms. 

For cold freeze-out evolutions [5] we use a second den- 
sity parameterization by Panov and Janka |14j . 

p(t) = p exp(-t/r) (2a) 

= ( 2b ) 

where po is the density at time t — and t$ is the switch 
over point corresponding to a temperature of Tg = 2. 
Equatio n |2a| gives early time behavior of the outflow and 
equation |2b| pertains to late time behavior. 

In modern r-process studies it is typical to approximate 
the thermodynamics by assuming the ejecta is radiation 
dominated. However, we require the full thermodynam- 
ics and so we instead use the solver from [51] to generate 
all trajectories. 

Our nucleosynthesis calculations consist of a multi- 
tiered algorithm which is coupled by three canonical 



stages. This nucleosynthesis network code was used in 
previous neutron capture studies by both Beun [2] and 
Surman [3J. 

The simulation starts with the regime of nuclear sta- 
tistical equilibrium. During this stage entropy and den- 
sity effectively determine all thermodynamic quantities. 
The second stage of the simulation employs an intermedi- 
ate reaction network [55J with PARDISO solver [551 157j . 
At this point the calculations include electromagnetic, 
strong and weak interactions. 

The third and final stage of our network calculation 
consists of a reduced r-process network as described in 
[481 158j . The primary reaction channels for nuclides in 
this section of the reaction network are beta-decay, neu- 
tron capture, and photo-dissociation. A detailed analysis 
regarding the formation of the rare earth elements during 
this stage is given in [5j [51] . 

In the third stage we track individual nuclear abun- 
dances by solving a set of differential equations given by 
the short hand notation (see [53] and pQ for details), 

Y(Z,A) = ]T \z>,a>Y z ,, a ,+ Y< pNA(<rv) z ,, A ,Y z ,, A ,Y n 

Z>,A> Z'.A' 

(3) 

where the quantity Y(Z,A) represents the time rate of 
change of abundance for nuclide (Z,A), Yz<, A < is the 
abundance for nuclide (Z',A'), p is the density, N A is 
Avagardo's number, (cry) is Maxwellian averaged neutron 
capture cross section for nuclide (Z' , A') and Y n the free 
neutron abundance. The first term on the right hand side 
holds information about beta decay modes and photodis- 
integrations. The rate, Xz>,A> could be one of following: 
the /3-decay rate with emission of j neutrons, \pj n or 
the photo-dissociation rate, A 7 . The second term on the 
right hand side includes reactions with neutrons. 

Nuclear flow is monitored by following individual terms 
in equation [3] For example, in the following equation we 
explicitly write out each term which gives the total flow 
in and out of a nucleus: 

Y in (Z,A) = \p(Z-l,A)Y(Z-l,A) 
+\pi n {Z -1,A + 1)Y(Z - 1, A + 1) 
+X P 2n(Z -1,A + 2)Y(Z - 1, A + 2) 
+\0 3n (Z -1,A + 3)Y(Z - 1, A + 3) 
+\ 1 (Z,A + l)Y(Z,A+l) 
+ (av) z , A -iY(Z,A-l)pN A Y n (4a) 

Y out (Z,A) = [X p (Z ) A) + Xp n (Z i A) 
+X, 32n (Z, A) + \p 3n (Z, A)]Y(Z, A) 
+A 7 (Z, A)Y(Z, A) + (av) z , A Y(Z, A)pN A Y n (4b) 

The difference of the two terms in equation[4]yields equa- 
tion [3l 

In general neutron reaction rates are of the following 



3 




120 130 140 150 160 170 180 190 200 
mass number 



FIG. 1: Shows the solar r-process abundance pattern, Nq <t versus atomic mass (data from [S5])- The three main regions of 
the pattern are highlighted. The A=130 peak is shaded white, the rare earth elements are shaded light gray and the A=195 
region is shaded in dark gray. Abundance scale is arbitrary. 



functional form (Fowler et al [60]): 




Ecr(E) exp{~E/k,T)dE 



(5) 

Neutron capture cross sections can vary by orders of 
magnitude between nuclear models [5] so we explore a 
variety of nuclear physics input: Finite Range Droplet 
Model (FRDM) [6THS3] . Extended Thomas- Fermi with 
Strutinsky Integral and Quenching (ETFSI) [62 -61] ver- 
sion 17 of the Hartree Fock Bogoliubov model (HFB17) 
[55] [70] and Duflo-Zuker (DZ) [5J]. The neutron capture 
rates of HFB17 and DZ were calculated using the TALYS 
code [ST]. The /3-decay rates used in our r-process net- 
work come from |68j . Calculations with different nuclear 
models show qualitatively similar behavior. So in our 
analysis we focus the discussion using one nuclear model, 
FRDM. We display results for all three nuclear models 
in section PVTl 

For the conditions we use entropy per baryon in units 
of Boltzmann's constant, S = 165, dynamical timescale, 
Tdyn = 85ms and electron fraction, Y e = .30 for a hot 
freeze-out and S = 105, Td yn — 50ms and Y e — .30 for a 
cold freeze-out. These conditions were chosen based on 
the procedure given in |69j . 

III. RESULTS 

A. Neutron Capture Rate Studies 

A neutron capture rate study consists of a 'baseline' 
simulation where the appropriate input of astrophysical 
and nuclear parameters are established. The final out- 
put is an abundance pattern (Y^ asehne ) as a function of 



atomic mass, A. Subsequent simulations are then con- 
ducted with the same input data but with the neutron 
capture rate of a single isotope in the rare earth peak 
changed by a factor of K = 5, 10, 50, or 100; each pro- 
ducing a final abundance pattern, Y K (A). 

A change in the neutron capture rate can lead to one 
of two distinct physical processes. A 'neutron capture ef- 
fect' occurs when change in capture rate results in path- 
way changes that shift material from the nucleus (Z, A) 
whose rate has been changed, to a nucleus (Z,A + 1), i.e. 
to the right on the NZ-plane. Or a 'photo-dissociation 
effect' which occurs when a change in capture rate on 
nucleus (Z, A) results in pathway changes that shift ma- 
terial from nucleus (Z,A + 1) to the nucleus (Z,A), i.e. 
to the left on the NZ-plane. We only find a neutron 
capture effect among the rare earth elements as photo- 
dissociation effects typically require highly populated nu- 
clei, see [3] for a detailed discussion of this effect. 

For each study, we consider the magnitude of influence 
a single neutron capture rate has on the final abundance 
pattern using the quantity F. 

Fk = ioo y lY {A) ~, A l - (6) 

A Y A 

A larger F value represents more deviation (percent 
change) from the baseline simulation. If F = then the 
neutron capture simulation abundance pattern and the 
baseline abundance pattern are equal; the neutron cap- 
ture rate change had no effect on the abundance pattern. 

We highlight the results of two different neutron cap- 
ture rate studies containing 126 nuclei in Figure [2] Each 
element's neutron capture rate has been changed by a fac- 
tor of K = 10 and the magnitude of its effect on the final 
abundances abundance is denoted by the degree of shad- 
ing. A darker gray (darker red online) represents a larger 
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effect. White boxes represent capture rates with little to 
no effect on the final abundance pattern (0 < F < 100). 
Each darker shaded nuclei represents an approximate fac- 
tor of two increase in the i^-measure with the darkest 
shade representing F > 400. Panel (a) of Figure [2] shows 
the effects of neutron capture rate changes under a hot 
evolution and panel (b) shows the effects of neutron cap- 
ture rate changes under a cold evolution. 



B. Comparing Hot and Cold Evolutions 

While the overall distribution of influential nuclei in 
Figure [2] is similar between the hot and cold freeze-out 
trajectories there are two prominent differences between 
these environments. First, there is a visible favoring of 
nuclei with even number of neutrons (even-N effect) in 
the cold evolution and an odd-N effect occurring at early 
times in the hot freeze-out evolution. Second, the mag- 
nitude of the neutron capture rate effect can vary for 
the same nuclei under the two r-process environments 
throughout the central shaded region in Figure [2j 

The r-process path is the time sequenced set of most 
abundant isotopes. It is useful to follow the path through 
the NZ-plane to understand the distribution of influen- 
tial nuclei across the rare earth region. In Figure [2j the 
path begins far from stability in the lower right corner of 
the NZ-plane progressively making its way back towards 
stability (top left) as the temperature falls and /3-decay 
begins to dominate the nuclear flow. 

As the path moves through the bottom right corner 
of the hot evolution (n, 7) T± (7, n) equilibrium is still 
in effect. Under these conditions small changes in neu- 
tron capture rates have no effect on the flow of material 
through the particular isotope. Due to equilibrium, the 
flow simply readjusts to compensate for the change. In 
the cold evolution, all photodisintegrations have frozen 
out and so the path is controlled by neutron captures 
and /3-decays only [H[S]. In this case, changes to neutron 
capture rates can impact final abundances. 

The path next encounters slower neutron capture rates 
in the central region of the figure. For those nuclei which 
are out of equilibrium and have significant abundance, 
changes in capture rate can now produce measurable ef- 
fects [5]. 

The nuclei in the top left portion are populated pri- 
marily via /3-decay. All other reactions, including neu- 
tron capture have frozen out. Thus we find little to no 
impact of neutron capture rates in this region. 



C. Changes to the Final Abundance Pattern 

We highlight the effect of an individual neutron cap- 
ture rate change has on the final abundances in Figure [3] 
These nuclides were chosen from the hot and cold neu- 
tron capture studies of Figure [2] Each neutron capture 



rate has been changed by a factor of K = 10 and the re- 
sultant final abundance pattern is compared to both the 
baseline and solar abundances. 

Inspection of Figure [3] shows that changes to neutron 
capture rates in the rare earth region produce only lo- 
cal changes in final abundances and that these changes 
are significant even for changes by a factor of K = 10. 
In both hot and cold evolutions 156 Nd shows a similar 
behavior under change in capture rate. This is not al- 
ways true for each nuclide, compare e.g. 168 Gd. We also 
find that individual nuclei exhaust their neutron capture 
effect at different values of K. For instance, in the hot en- 
vironment, the 165 Eu neutron capture effect is exhausted 
near K = 50 while in the cold environment, the 165 Pm 
neutron capture effect is maximal near K = 10. 

IV. ANALYSIS OF THE NEUTRON CAPTURE 
EFFECT 

We now consider the factors which influence the mag- 
nitude of the neutron capture effect. A necessary condi- 
tion for large neutron capture effect is that the nuclide 
exhibit large nuclear flow through the neutron capture 
channel. In order for a large flow through the neutron 
capture channel to have a significant influence on the 
abundance pattern the neutron capture flow must also 
be unsaturated. Flow saturation, which we denote as S, 
occurs when the sum of material flowing through input 
channels matches the flow of a single output channel. In 
the case where the output channel is neutron capture 
the channel is 'saturated' (S — 0) when it is matched by 
photo-dissociation and /3-decay in-flows. We use these 
concepts to explain the relative differences seen in F be- 
tween nuclides in the hot and cold evolutions. 



A. Large Nuclear Flow 

A large neutron capture flow means significant mate- 
rial transportation through the neutron capture reaction 
channel. The last summand in equation |4b| contains the 
relevant information on the movement of material via 
neutron capture flow out from isotope {Z, A) to isotope 
(Z.A + 1). We provide it for convenience: 

neutron capture flow = (av) z ^Y(Z, A)pNaY„ (7) 

This equation consists of three main ingredients: 
the thermally averaged neutron capture cross section, 
(av) z,A, the abundance of the particular isotope, Y(Z, A) 
and the neutron number density, pN aY u . The interplay 
between these three components determines the size of 
the neutron capture flow. 

For the energy ranges explored in the two astrophysi- 
cal environments considered here, {<jv)z,a is a relatively 
flat function of temperature. Therefore, for a fixed nu- 
clear data set the differences between classical and cold 
neutron capture flows is primarily due to differences in 
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FIG. 2: Shows neutron capture rates that significantly influence the abundance pattern. The results of two separate neutron 
capture rate studies are shown for a hot freeze-out, panel (a) and cold freeze-out, panel (b). In both cases individual neutron 
capture rates were changed by K = 10. Darker shades represent an increased effect on the abundance pattern. In order of 
lightest to darkest each shade represents: white (0 < F < 100), light (100 < -F < 200), medium (200 < F < 400), darkest 
(F > 400). Above the dotted line neutron capture flows are not large enough for an increase in a neutron capture rate to 
produce a significant neutron capture effect. 
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FIG. 3: Shows the effect of particular neutron capture rates on the rare earth abundances. Simulations were performed using 
hot (panel a) cold (panel b) freeze-out evolutions. The baseline curve, Y basel ' ne is represented by a bold black line and the 
solar data is repressented by a solid gray line. For both types of trajectories we show five curves, Y K=1 ° each representing a 
simulation where a single neutron capture rate was changed by a factor of 10. 
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FIG. 4: Compares neutron capture flows in baseline simulations for the hot (solid) and cold (dashed) evolutions as a function 
of simulation time. Also shown are the neutron capture rate, nuclear abundance and neutron number density. The product 
of these three determines the neutron capture flow. The neutron capture rate change has the most influence around the time 
when 5* is maximal (star). The left panel highlights 168 Gd and the right panel 171 Dy. 



the neutron number density and elemental abundance. 
Both components are significantly influenced by the as- 
trophysical environment and nuclear data set. Figure 
[4] shows the interplay between the three components of 
equation [7] as a function of simulation time for the base- 
line case. Two nuclides are highlighted, 168 Gd in the left 
panel and 171 Dy in the right panel. 

To maximize the neutron capture rate effect we search 
for large out of equilibrium neutron capture flows in the 
baseline simulation. In order to measure the magnitude 
of the neutron capture flow in the baseline simulation we 



compute, 

r _ J (<rv)z,AY(Z,A)pN A Y n dt 

L ~ Y*(Z,A) (8) 

where Y*(Z,A) is the abundance when the nucleus is 
furthest from saturation (5 maximal), the integrand is 
equation [7] and the integral is taken over simulation time. 

When L > .2, flow through the neutron capture chan- 
nel is sufficiently large for a change in neutron capture 
rate to produce significant change in the final abun- 
dances. We highlight this in both panels of Figure [2] by a 
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dotted line. Above the dotted line, L < .2 and below the 
dotted line L > .2. The measure L shows quantitatively 
how changes in the neutron capture rates of nuclei above 
the dotted line have no influence on the final abundances. 



B. Flow Saturation 

In the region of Figure [2] where L > .2, we can use flow 
saturation to understand the differences in the magnitude 
of the neutron capture effect among nuclei. 

To measure saturation in the neutron capture channel 
we take equation |4a| and subtract equation [7J 

S = Y in (Z, A) - (av) ZtA Y(Z, A) P N A Y n (9) 

Saturation occurs when 5 = 0. Changes to a neutron 
capture rate under saturation have no effect on the flow 
of material because the output channel is limited by the 
in- flowing channels. When S < 0, more material is flow- 
ing out through the neutron capture channel than is flow- 
ing into the input channels. When S > 0, the neutron 
capture flow is smaller than the in-flowing channels. 

For a large neutron capture effect it is crucial that S is 
large and greater than zero so that the neutron capture 
flow is furthest from saturation. We can approximate 
the time at which a large neutron capture flow is impor- 
tant in producing a large neutron capture effect in each 
simulation by finding the maximum of S. 

When predicting the magnitude of the neutron capture 
effect it is useful to define the integral of equation [9] over 
simulation time for S > 0, 

S= J Sdt (10) 

s>o 

which we call the unsaturated index. Physically, it is the 
amount of material that is flowing into the nucleus but 
is not flowing out via neutron capture. Thus, this ma- 
terial could be directed out through the neutron capture 
channel if the neutron capture rate were increased. 

A special case of neutron capture flow saturation oc- 
curs as nuclei fall out of (n, 7) ^ (7, n) equilibrium 
in hot scenarios. If the temperature is high enough 
photo-dissociation rates are large so that the net flow 
to the right can be limited by leftward flowing material 
in the photo-dissociation channel. If the nuclide is in 
(n, 7) ^ (7, n) equilibrium then the photo-dissociation 
and neutron capture terms in equation [9] are large and 
cancel so that S sa 0, see bottom right corner of panel 
(a) of Figure [2} 

Odd-N nuclei are particularly susceptible to flow satu- 
ration as they have smaller separation energies than even- 
N nuclei. This means that for odd-N nuclei the neutron 
capture photo-dissociation rate pair (Z, A) and (Z, A + l) 
tends to fall out of equilibrium sooner than the rate pair 
(Z, A— I ) and (Z, A) [3]. Thus, the F-measure is sensitive 
to neutron capture rates on odd-N nuclei far from stabil- 
ity in a hot freeze-out. For example, in the top panel of 



Figure § 161 Ce, 165 Nd, 166 Pm and 170 Pm all fall out of 
equilibrium earlier than the surrounding nuclei and each 
exhibit a small neutron capture effect limited by flow sat- 
uration. 

Flow saturation also occurs in cold environments. In 
this case, neutron capture rate effects are limited for odd- 
N rare earth isotopes because they have faster neutron 
capture rates and /3-decay rates compared to even-N rare 
earth isotopes. Faster neutron capture rates in odd-N nu- 
clei means the first term in equation[9]is larger for even-N 
nuclei than for odd-N nuclei. Faster beta decay rates in 
odd-N nuclei also implies the first term in equation [9] is 
larger for even-N nuclei than for odd-N nuclei. The net 
effect is that neutron capture flows of odd-N nuclei are 
closer to saturation (S « 0) than the flows of even-N 
nuclei. Therefore, even-N nuclei are favored by the F- 
measure, see bottom right corner of panel (b) of Figure 
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C. Flow Saturation As A Predictor For The 
Magnitude Of The Neutron Capture Effect 

Computing the unsaturated index using the baseline 
simulations, we can now predict and understand the rel- 
ative differences observed in the F-measure of the same 
nuclei in the central region of Figure [2] between hot and 
cold freeze-out conditions. 

For 168 Gd the unsaturated index is larger in the hot 
baseline, S — 6.73 x 10 -3 than in the cold baseline, 
S = 1.27 x I0~ 3 . This implies the neutron capture effect 
should be larger in the hot scenario as can be verified by 
comparing the two panels in Figure [2j Saturation also 
estimates when the neutron capture effect is important, 
S maximal. Returning to Figure [ZJ we see that this point 
(star on the figure) is sensitive to thermodynamic condi- 
tions and may occur when the neutron capture flow is 
not at its largest value. 

The neutron capture flow components of Dy near 
the time of maximal S also leads to a large neutron cap- 
ture flow, albeit a magnitude smaller than the neutron 
capture flow of 168 Gd. For m Dy the unsaturated index 
is larger in the cold baseline, S = 2.44 x I0~ 4 as com- 
pared to the hot baseline, S — 1.23 x I0~ 4 . This implies 
the neutron capture effect should be larger in the cold 
scnairo and can again be verified by comparing the two 
panels in Figure [2j 

Explaining the relative differences seen in F by com- 
paring unsaturated indices and then confirming the result 
with F-measure values works extremely well. However, a 
word of caution is necessary: The F defined in equation[6] 
is a sum of percent abundance differences and the unsatu- 
rated index, equation |10| contains total abundance yield 
information. Thus if one prefers exact agreement be- 
tween the result of the capture study (F-measure value) 
and the unsaturated index, one should only use differ- 
ences in abundance or massfraction as the F-measure 
(for example, see F defined in [3]). 



V. ANALYSIS: PATHWAY CHANGES 

In the previous sections we have discussed the magni- 
tude of the neutron capture rate effect and demonstrated 
how changes in neutron capture rates can influence the 
final abundance distribution. We now turn our attention 
to the pathway changes produced by the rate changes, 
focusing on individual nuclei in the K = 10 case. 

In order to study changes in the path, for each element 
in the network we produce a real number AY(Z, A) rep- 
resenting the total change of the abundance yield over 
r-process simulation time, 

AY(Z, A) = J [Y£J(t) - Y h z a X Une [t)\dt (11) 

where Y^ c ^ is the abundance of isotope (Z, A) when one 
neutron capture rate has been changed and yj a ^| e ' me i s 
the abundance of isotope (Z, A) in the baseline simu- 
lation. If AY(Z, A) is positive then more material re- 
sides in the particular nucleus during the capture rate 
simulation than during the baseline simulation. Con- 
versely, less material resides in the nucleus during the 
capture rate simulation than during the baseline simu- 
lation if the value of AY(Z,A) is negative. For capture 
rate changes in the rare earth region most nuclei in the 
network have AY(Z,A) = 0, except for locally around 
the nucleus whose capture rate has been changed. 

For each neutron capture rate study we normalize the 
set of Ai^'s and represent their magnitude by colors in 
the NZ-plane of Figures [5] and [6j Darkest solid shades 
(online dark green) represent the largest order of magni- 
tude positive change while light solid shades (online light 
green) represent the next largest positive change. Dark- 
est hatches (online dark red) represent the largest order 
of magnitude negative change while light hatches (online 
light red) represent the next largest negative change. 

Pathway changes can vary between the same nuclei 
under different simluations for a number of reasons. This 
includes variations due to astrophysical conditions, the 
onset of r-process freeze-out, availability of free neutrons, 
large neutron capture flow or flow saturation. 

Of the ten elements shown in Figure [3] we select four 
of them for a study of the pathway changes. We examine 
Neodymium-156 in the context of pathway changes due 
to differences in astrophysical conditions. We then con- 
sider pathway changes of two different nuclei Europium- 
165 and Promethium-165 to highlight the importance of 
large neutron capture flow of surrounding nuclei. 

In Figure[5]we see the changes in nucleosynthetic path- 
ways of Neodymium-156, which are slightly greater in the 
classical senario than in the cold senario. In both cases, 
before changing the neutron capture rate 156 Nd was pop- 
ulated from a /3-decay channel from 156 Pr. Because of 
the quick /3-decay rate of 156 Nd the material continued 
to flow to 156 Pm; following the dotted arrows. 

After increasing the neutron capture rate by a factor 
of K = 10, under the hot trajectory (left panel), the neu- 
tron capture channel is enhanced. The flow of material 



now travels through the neutron capture channel rather 
than the /3-channel resulting in material being deposited 
in new elements as highlighted in a light solid shade; fol- 
low solid arrows. 

Under the cold evolution we find a similar effect. How- 
ever, the neutron capture rate of 156 Nd in the cold en- 
vironment baseline is half the value of the capture rate 
of 156 Nd in the hot environment baseline. At the time 
when S is maximal, less material will travel through the 
capture channel in the cold senario with increase in the 
capture rate resulting in a slightly more constrained path- 
way. We see this in Figure [5] by comparing panel (a) 
and (b). In both environments 156 Nd reaches saturation 
around if = 10. 

In Figure [6] the nucleosynthetic pathway changes of 
Europium-165 are displayed for the hot evolution (left 
panel) and the pathway changes of Promethium-165 for 
the cold evolution (right panel) under a capture rate 
change of K = 10. These two nuclei's capture rates 
have similar effects on the abundance pattern, F « 400, 
but the pathway changes for 165 Pm extend through nine 
units of atomic mass while the pathway changes of 165 Eu 
extend through only four. This discrepancy arises from 
the differences in the flow through the neutron capture 
channel of the surrounding nuclei. 

In the case of 165 Eu the capture rate change is impor- 
tant during relatively late times. At this point in the hot 
simulation the temperature has fallen drastically and free 
neutrons are relatively scarce. The unchanged neutron 
capture channels of nuclei surrouding 165 Eu have diffi- 
culty competing with the increasing /3-decay flow in the 
region. The new pathway is thus limited from branching 
out resulting in most material being shifted to 166 Gd. 

For 165 Pm the capture rate change is important farther 
from stability (at earlier times) in the cold evolution and 
it also does not have to compete with photo-dissociation 
flows. Note that 165 Pm has no effect on the abundance 
pattern under the hot evolution, see bottom panel of Fig- 
ure [2] In addition, 165 Pm is an even-N nucleus so that 
it is far from saturation compared to odd-N nuclei popu- 
lated at early times. At this point in the cold simulation 
the unchanged neutron capture rates of the surrounding 
nuclei can still compete with /3-decay rates in the region 
resulting in the formation of many secondary pathways 
beyond 166 Pm. 



VI. NEUTRON CAPTURE RATE DECREASES 
& NUCLEAR PHYSICS UNCERTAINTIES 

For completeness, we also conducted neutron capture 
rate studies where the rates were decreased by factors of 
5, 10, 50 and 100. The same general analysis presented 
above can be applied to neutron capture rate decreases. 
The major difference arises when saturated nuclei (in the 
baseline simulation) are moved further from saturation 
by decreasing neutron capture rate, reducing the second 
term in equation [9j 
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FIG. 5: Shows the change in the nucleosynthetic pathway when the neutron capture rate of 156 Nd is changed. In the baseline 
simulation the flow of material through 156 Nd occurs primarily in the /3-decay channels. When the neutron capture rate is 
increased by a factor of K — 10, the flow out of 156 Nd is primarily through neutron capture. This occurs in both types of 
trajectories hot (left panel) and cold (right panel). Jagged arrows represent nuclear flow in the baseline simulation while solid 
bold arrows represent flow with the changed capture rate. Relative decreases in abundances along the path are denoted by 
hatched gray tones (online red) while relative increases are represented by solid gray tones (online green). 
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FIG. 6: Shows the change in nucleosynthetic pathway when the neutron capture rate of 165 Eu (left panel) and 165 Pm (right 
panel) are modified by a factor of K = 10. In the left panel we study the 165 Eu neutron capture effect under the hot trajectory. 
In the right panel we study the 165 Pm neutron capture effect under the cold trajectory. In both cases the baseline simulations 
exhibit flow dominated by /3-decay. The neutron capture effect of 165 Eu changes the pathway resulting in most material being 
deposited in 166 Eu and 166 Gd. The neutron capture effect of 165 Pm represents an extreme case where multiple new channels 
are opened. The markings are the same as in Figure [5] 



To summarize the effects of nuclei whose neutron cap- 
ture rates can significantly impact abundance patterns 
we generated a combined data set which includes four 
nuclear models (FRDM,ETFSI,HFB17,DZ) and two as- 
trophysical conditions (hot, cold) studying both increases 
and decreases by factors of 5, 10, 50 and 100. We sepa- 
rated the data based on whether the rate was increased 
or decreased. We then selected the maximal F for each 
nuclei across these studies, the results of which are shown 
in Figure [7j 

Panel (a) shows those nuclei whose neutron capture 
rates can significantly impact final abundances when the 



rates are increased. The lengthy freeze-out phase of cold 
evolutions contributes to most of the shaded nuclei in the 
bottom right portion of the panel. Here, Even-N nuclei 
are favored across different nuclear data sets due to flow 
saturation in odd-N nuclei. The neutron capture effect of 
nuclei near the rare earth peak, with = 102, N = 104 
and = 106 and from Praseodymium up to Europium 
is exhausted for changes in neutron capture rates of only 
K = 5. Capture rate changes above K = 5 do not matter 
for these nuclei as their neutron capture channel becomes 
exhausted due to saturation. 

Panel (b) shows those nuclei whose neutron capture 
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rates can significantly impact final abundances when the 
rates are decreased. We find that decreases in neutron 
capture rates in the rare earth region tend to have smaller 
effects than increases. Because neutron capture rate de- 
creases restrict the flow through the neutron capture 
channel we find that larger neutron flow, L > .8, is re- 
quired to produce significant change to the final abun- 
dance pattern. Thus, the distribution of influential neu- 
tron capture rates is shift to more neutron rich nuclei. 

VII. CONCLUSIONS 

We have demonstrated the importance of understand- 
ing individual neutron capture rates in the rare earth re- 
gion of the r-process abundance pattern and shown that 
the distribution of influential nuclei can be elucidated by 
the concepts of large nuclear flow and flow saturation. 
These concepts are applicable across a variety of astro- 
physical conditions and nuclear models. An influential 
neutron capture rate leads to a 'neutron capture effect' 
where a change is effected in the abundance of nearby 
higher A nuclei. 

Many nuclei show significant leverage on the final 
abundances with small neutron capture rate change (by 
a factor of K — 5) , as shown in Figure [7j The distri- 
bution of important rare earth neutron capture rates in 
the NZ-plane occurs in a narrow band approximately 10 
to 20 neutrons from stability. The overall distribution 
is remarkably similar across a variety of freeze-out con- 
ditions and differing input nuclear physics, e.g Figure [7] 
and is in agreement with the prediction of the location of 
important neutron capture rates in the rare earth region 
based on formation arguments [5]. 

To understand the magnitude of the neutron capture 
effect we introduced two concepts: (1) large nuclear flow 
and (2) flow saturation. Large nuclear flow means signifi- 
cant material transportation in a given reaction channel. 
The requirement of large nuclear flow in the neutron cap- 
ture channel tends to rule out nuclei that are less than 10 
neutrons from stability in the rare earth region as these 
nuclei are populated primarily via /3-decay. Flow satu- 
ration occurs when the sum of material flowing through 
input channels matches the flow of a single output chan- 
nel. In the case where the output channel is neutron 



capture, the channel is 'saturated' when it is matched by 
photo-dissociation and /3-decay in-flows. Changes to a 
neutron capture rate under saturation have no effect on 
the flow of material because the output channel is lim- 
ited by the in-flowing channels. In like manner, changes 
in capture rates by very large factors are also forbidden 
because the input channels become exhausted. 

Flow saturation is useful in understanding the details 
of the pattern of influential nuclei. For instance, un- 
der hot freeze-out conditions, photo-dissociation rates are 
large so that the net flow to the right in the NZ-plane 
can be limited by leftward flowing material in the photo- 
dissociation channel. This limits the neutron capture ef- 
fect for nuclei far from stability cither in (n, 7) <=^ (7, n) 
equilibrium or for those nuclei just coming out of equi- 
librium because the channel is saturated. Under cold 
freeze-out conditions, odd-N nuclei are closest to satura- 
tion since they tend to have faster neutron capture rates 
and /3-decay rates. Thus, large neutron capture effects 
are generally found in even-N nuclei. 

The concepts of large nuclear flow and flow satura- 
tion are general concepts which are applicable beyond the 
scope of individual neutron capture studies. We expect 
that the overall distribution of influential nuclei would 
still exhibit the same qualitative behavior, if for instance, 
one chose to study groups of neutron capture rates. 

In order to disentangle the information contained in 
the abundances found in nature, detailed knowledge of 
nuclear physics including masses, /3-decay rates and neu- 
tron capture rates must be known. The concepts of large 
nuclear flow and flow saturation introduced here are es- 
sential concepts for studying any reaction rates relevant 
to the r-process. 
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FIG. 7: Identifies all nuclei in a combined data set whose neutron capture rates can significantly impact the final abundance 
pattern when rates are (a) increased and (b) decreased. The data set includes calculations using a combination of nuclear 
models (FRDM,ETFSI,HFB17,DZ) anc l conditions (hot and cold). Shaded nuclei have F max ~ 200 or more abundance change 
for a neutron capture rate increase or decrease of K = 5 (darkest shaded squares), K — 10 (medium shaded squares), K — 50 
(light shaded squares) and K = 100 (lightest shaded squares). For each nuclei the maximal F was chosen among the data sets. 
Nuclei shaded white never produce a significant effect (F > 200) under any data set for any rate change, K. 
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